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Temperature  Using  Embedded  Temperature  Sensors 

(Preprint) 

Edward  B.  Coy1 

Air  Force  Research  Laboratory,  Edwards  AFB,  CA,  93524 

Approximate  solutions  have  been  obtained  for  surface  temperature  and  heat  flux  for  the  case  of 
two  embedded  temperature  sensors.  The  solutions  have  been  verified  and  the  range  of  validity  has 
been  established  using  several  methods  including  comparisons  with  an  exact  analytical  solution 
for  a  linear  problem  and  a  numerical  calculation  for  a  non-linear  problem.  The  performance  of  the 
model  is  presented  in  both  the  frequency  domain  and  the  time  domain.  A  propagation  of  error 
analysis  is  presented  and  is  used  to  establish  the  optimum  spacing  between  the  sensors.  The 
solutions  place  no  restrictions  on  the  boundary  or  initial  conditions  and  rely  only  on  current  values 
of  temperature  and  its  rate  of  change.  The  temperature  dependence  of  transport  properties  is 
accounted  for  in  an  approximate  way.  The  method  requires  minimal  computation  and  is  suitable 
for  implementation  in  a  real-time  sensor. 


Nomenclature 


Fo 

k 

cp 

Tt 

f, 

t 


=  Fourier  number,  atjx 1 
=  Thermal  conductivity,  Watt/meter2-Kelvin 
=  Fleat  capacity,  Joule/Kilogram-Kelvin 
=  Temperature,  Kelvin 

=  Rate  of  change  of  temperature,  Kelvin/second 
=  Time,  second 
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=  Heat  Flux,  Watt/meter2 
=  Distance  from  heat  flux  surface,  meter 
=  Thermal  diffusivity,  meter2/second 
=  Phase  Angle,  radian 
=  Angular  frequency,  radian/second 


I.  Introduction 

The  measurement  of  transient  heat  flux  and  surface  temperature  in  heat-sink  combustion  chambers  continues  to 
present  technical  challenges  to  the  instrumentation  engineer.  Sensor  failure  rates  are  high  and  measurement 
accuracies  and  uncertainties  are  not  well  characterized.  These  shortcomings  have  had  a  significant  impact  on  some 
recent  programs  which  have  used  heat  sink  test  articles  to  acquire  data  for  the  validation  of  heat  transfer  predictions 
at  liquid  rocket  engine  operating  conditions1'2. 

There  are  numerous  types  of  heat  flux  sensors  but  a  relatively  small  subset  is  capable  of  operating  in  rocket 
chamber  conditions  where  heat  flux  levels  can  exceed  108  W/m2  and  surface  temperatures  of  1000  K  are  typical. 
Diller3  reviewed  the  devices  that  have  been  used  and  organized  them  into  methods  that  rely  on  temperature 
differences  over  a  spatial  distance  with  known  thermal  resistance  and  temperature  differences  over  time  with  known 
thermal  capacitance.  The  most  commonly  used  method  has  been  the  coaxial  thermocouple,  which  is  an  example  of 
the  second  type.  A  thermocouple  junction  is  formed  on  the  surface  of  the  chamber  between  a  wire  of  one  type  of 
thermocouple  material  and  a  surrounding  sheath  of  another  type.  The  heat  flux  is  determined  from  the  measured 
temperature  boundary  condition  using  a  one -dimensional  transient  solution  to  the  heat  equation.  The  junction  is 
typically  very  thin  and  is  often  formed  by  lightly  scratching  the  surface  to  drag  filaments  of  one  type  of  material 
across  the  electrically  insulating  layer  to  the  other  type.  In  some  applications,  when  erosion  of  the  surface  occurs, 
the  junction  is  continuously  reformed  and  this  has  led  to  the  description  of  coaxial  thermocouples  as  “eroding 
thermocouples.”  However,  in  heat-sink  chambers  it  is  quite  common  to  find  that  the  junction  disappears  at  some 
point  during  a  test  and  the  sensor  fails. 

Other  methods  have  been  developed  which  do  not  rely  on  surface  temperature  measurements  but  embed  the 
sensors  within  the  wall  where  they  are  protected  from  erosion.  In  the  null-point  calorimeter  a  hole  is  drilled  from 
the  backside  of  the  chamber  wall  and  a  thermocouple  is  inserted.  The  bead  is  brazed  or  resistance  welded  to  the 


2 

American  Institute  of  Aeronautics  and  Astronautics 


bottom  of  the  hole.  “Null-point”  refers  to  a  distance  from  the  bottom  of  the  hole  to  the  inner  wall  where  the 
disturbance  to  the  flow  of  heat  caused  by  the  hole  results  in  the  junction  reading  nearly  equal  to  the  inner  wall 
temperature.  The  construction  of  null-point  calorimeters  is  challenging.  The  junction  cannot  be  visually  inspected 
and  large  measurement  errors  can  result  from  manufacturing  flaws4. 

Another  method  using  embedded  temperature  sensors  is  the  plug-type  heat  flux  gauge  of  Liebert5.  An  annular 
groove  is  machined  into  the  chamber  wall  to  form  a  post  and  thermocouples  are  attached  at  several  axial  points 
along  the  outside  of  the  post.  A  polynomial  curve  is  used  to  extrapolate  the  temperatures  to  the  wall  position  and  an 
integral  method  is  used  to  calculate  the  total  heat  load  to  the  plug  from  transient  temperature  measurements.  Two- 
dimensional  effects  can  be  significant  in  this  type  of  device.  The  dimensions  of  the  groove  are  critical  and 
significant  errors  can  result  from  the  disturbance  to  the  flow  of  heat6. 

Recently,  Conley,  Vaidyanthan  and  Segal  used  embedded  temperature  sensors  to  measure  transient  heat  flux  in  a 
subscale  combustion  chamber7.  Transient  temperature  data  were  reduced  to  heat  flux  using  an  equation  developed 
from  an  energy  balance  on  a  control-volume  between  the  sensors.  The  equation  included  terms  for  Fourier’s  law  of 
heat  conduction  and  heat  storage  and  nominal  constant  values  were  used  for  the  thermal  transport  properties.  The 
accuracy  of  the  method  was  estimated  to  be  +/-  2%  for  the  conditions  of  the  study;  however,  a  general  discussion  of 
the  accuracy  and  limitations  of  the  method  was  not  provided. 

In  this  paper  a  method  for  calculating  the  surface  boundary  conditions  in  transient,  heat-sink  test  articles  using 
embedded  temperature  sensors  is  described.  The  method  has  been  used  in  our  laboratory  over  the  past  two  years 
and  has  provided  satisfactory  solutions  to  the  major  technical  challenges  of  accuracy,  uncertainty,  time  response, 
failure  rate,  and  producibility.  The  one-dimensional  planar  case  is  presented  but  the  technique  can  be  also  be 
formulated  for  cylindrical  or  spherical  geometry. 

II.  Derivation  of  Model 

In  the  following  analysis  we  adopt  the  approach  of  Liebert  of  extrapolating  a  polynomial  to  the  surface.  The 
temperature  profile  in  the  body  (see  figure  1 .)  is  approximated  using  a  power  series  in  x  with  time  dependent 
coefficients. 

T(x,t)  =  YjCi(t)x"  (1) 

i= 0 
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However,  in  addition  to  matching  the  temperatures  at  the  measurement  points,  the  polynomial  is  also  required  to 
satisfy  the  heat  equation.  Specifically,  the  polynomial  is  required  to  match  the  second  derivatives  of  temperature 
with  respect  to  position  at  the  measurement  points.  An  estimate  that  is  adequate  under  most  circumstances  can  be 
obtained  by  neglecting  gradients  in  the  thermal  conductivity  and  using  the  linear  form  of  the  heat  equation. 

d2T  1  dT 

-^Y  =  ~^7  (2) 

ox  a  ot 

The  C  j  are  obtained  by  substituting  (1)  into  the  left  side  of  (2)  and  then  solving  for  the  C;  in  terms  of  the  measured 

temperatures  and  their  rates  of  change.  For  two  measurement  points  there  will  be  four  measured  quantities  and  the 
power  series  can  be  carried  to  the  fourth  term.  If  second  derivatives  are  available,  the  series  can  be  carried  to  the 
sixth  term  and  so  on. 


Evaluating  the  series  at  the  surface,  x=0,  and  combining  terms  in  order  to  make  the  linear  dependence  on  the 
measurements  apparent,  the  temperature  and  heat  flux  expressions  can  be  written  as  follows. 
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(3) 


T  =T 

1o  1i 


+  Tn 


d2Tx  xxx2(2x2  -jCj)  d%  xxx2(x2  -  2xx ) 


x-,  -  x, 


x,  -  X, 


cbr 


6  (x2  -  xx ) 


obr 


6  (x2  -  xx ) 


<7o  =~k(T0) 


Tt 


-1 


■  +  T ' 


d  Tx  x[  -2xxx2  -2x2  +d~T2  2xx  +2xxx2  -x2 


2  \ 


dxl 


6(x2  -  xx) 


dx2 


6{x2  -  x, ) 


(4) 


In  deriving  these  expressions,  no  assumptions  have  been  made  with  respect  to  the  boundary  condition  at  the  top 
of  the  block  and  in  fact  it  is  completely  arbitrary.  For  example,  the  top  surface  could  be  actively  cooled.  Also  the 
second  sensor  at  x2  does  not  need  to  be  embedded  within  the  block  and  can  be  located  at  the  top  surface  and  the 
sensor  at  X/  could  be  located  at  the  x=0  surface.  Furthermore,  we  have  made  no  assumptions  with  respect  to  the 
origin  of  time;  the  solution  relies  only  on  the  current  values  of  temperature  and  rates  of  change  of  temperature. 
Referring  back  to  the  categorizations  of  methods  defined  by  Diller  ',  this  method  is  a  hybrid  of  the  two  types  as  it 
relies  on  both  spatial  and  temporal  variations  in  temperature. 

The  method  used  to  evaluate  the  rate  of  change  of  temperature  is  critical  to  the  success  of  the  technique.  A 
simple  finite  difference  calculation  will  most  likely  result  in  an  unacceptable  level  of  noise  in  the  derivative.  A 
polynomial  smoothing  filter,  also  known  as  a  Savitzky-Golay8,9  filter,  can  be  effective  in  reducing  noise  and  can  be 
calculated  efficiently.  The  filter  is  implemented  as  a  convolution  in  the  time  domain. 


Ti(fn)=  YjajTi^n+j) 
j=-m 


(5) 


r,(M=  YJb,nt.„)  (6) 

j=-m 


Substituting  (5)  and  (6)  into  (2),  (3)  and  (4)  we  arrive  at  the  final  forms  for  the  surface  temperature  and  heat  flux 
expressions. 
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(7) 


m 

T0(tn)=  ' Y.eiT,(tn+i)  +  fiT2(tn+i ) 

i=—m 


m 

q0  (tn  )  =  -k(T0 )  ^  g,Tt  (tn+l )  +  h,T2  (tn+l ) 

i=—m 


(8) 


The  coefficients  in  (7)  and  (8)  are  defined  as  follows. 


x~ 


Jf,  —  X, 


xlx2  (2x2  —x1) 
6a j  (x2  —  xl ) 


(9) 


fi  =  a. 


*i  1  b  X1x2(x2  —  2xt) 
x2  —  xl  '  6«2  (x2  —  xy ) 


(10) 


Si  =  a, 


1  ,7-^i  2x2 


x2  —  Xj  '  6a1(x2  —  Xj) 


(11) 


,  1  ,  2x,2  +  2x,  x2  —  x2 

hj  =  at - h  b:  1  1  2  2 


X- 


—  x,  !  6«r2  (x2  —  x, ) 


(12) 


Note  that  once  the  coefficients  have  been  calculated,  only  4m+2  multiplications  and  2m+l  additions  are  required 
to  evaluate  T0  and  q0 ,  or  about  100  processor  operations  are  required  for  a  typical  value  for  m  of  10. 

The  uncertainties  ST{)  and  Sq0  can  be  evaluated  by  treating  each  temperature  measurement  as  an  independent 

random  variable  and  using  the  “square  root  of  the  sum  of  squares”  approach.  In  the  following  we  assume  that 
uncertainties  in  thermocouple  locations  and  material  properties  are  negligible. 


«i«„) 


I 


\  1/2 


e,.2  ST2 


+  fi  Hi 


\i=—m 


j 


(13) 
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/  \  1/2 

'  m  \ 

aj0(tn)=k(T0)\  '£ig.2srl2+hi2sr22 

\i=—m  J 


(14) 


Minimum  values  of  ST(]  and  Sq0  can  be  achieved  by  making  an  optimum  choice  for  the  order  of  the  polynomial 


used  to  obtain  the  a ,  and  bl .  For  a  fixed  level  of  uncertainty  in  the  measurements,  the  coefficients  of  a  high  order 
method  will  result  in  a  larger  value  for  the  total  uncertainty;  however,  in  the  presence  of  high  frequency  content,  a 
high  order  polynomial  will  fit  the  trend  more  closely  and  reduce  the  values  of  8 Tn  ,  thereby  reducing  the  total 

uncertainty.  The  optimum  order  can  be  found  by  increasing  the  order  until  a  minimum  in  the  uncertainty  is  found. 
However,  this  increases  the  amount  of  computation  and  in  practice  we  find  that  the  quadratic  smoothing  filter  is  a 
good  compromise  and  can  be  used  as  the  default  method.  The  coefficients  for  the  quadratic  smoothing  filter  are  the 
following. 


Cl:  = 


£j2  2 ]jA 

j=—m  j——m 

m 

I  f 

j=-m 


2m  + 1 


If  I r 

j=—m  j=—m 


(15) 


m 

AfX./2 


In  applications  where  the  linearized  heat  equation  (2)  is  not  sufficiently  accurate,  improved  estimates  for  the 
C  j  coefficients  in  (1)  that  account  for  gradients  in  thermal  conductivity  can  be  obtained  by  a  second  iteration.  The 
heat  equation  for  this  case  can  be  written  as  follows. 


d2T  1  dT  1  dk 

' dT 

\ 

dx 2  a  dt  k  dT 

v  dx 

linear  J 

(17) 
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In  (17),  dT/dx  is  estimated  using  the  derivative  of  (1)  based  on  the  initial  estimates  for  theC;,  and  dk/dT  is 
obtained  from  a  correlation  for  the  material.  The  process  of  iteratively  improving  (17)  can  be  repeated  until  a 
convergence  criterion  is  satisfied.  The  new  estimates  for  d~T j ' dx~  can  then  be  used  in  (3)  and  (4)  directly. 


III.  Verification  of  Model 

To  illustrate  the  behavior  of  the  approximate  solutions,  we  make  comparisons  with  the  exact  solution  for  the 
temperature  within  a  semi-infinite  slab  exposed  to  an  oscillating  heat  flux  at  the  surface. 


qQ  =  q  cos  (cot) 


(18) 


T(x,t)  =  T0  +  y 
k 


—  /  \l/2  f 

q '  "  ' 


a 

\co) 


exp 


/  \l/2\  C 

'  CO  1 

cosl 


v  y 


cot  -  A| 


C  \ln  ^ 

'co  }  n 

4 


\2a  j 


(19) 


When  (19)  and  its  time  derivative  are  substituted  into  (2),  (3)  and  (4)  the  results  are  linear  combinations  of  sine 
and  cosine  terms  of  the  form: 


T0  =  Acos(cot  +  cj)x )  +  Bcos(cot  +  <f>2)  +  C  sin(&>f  +  cj>x)  +  Dsm(cot  +  <f>2) 


(20) 


The  four  terms  can  be  combined  into  a  single  term  of  the  following  form. 


T,  =  Gain  — 
k 


-  c  \1/2 
Vo 


a 

\<oj 


sin(<z>f  +  Phase ) 


(21) 


The  gain  and  phase  behaviors  are  shown  in  figures  2  and  3.  For  greater  generality,  the  results  are  plotted  as 
functions  of  two  non-dimensional  numbers.  The  abscissa  variable  is  a  dimensionless  angular  frequency  which  is  the 
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reciprocal  of  a  Fourier  number  based  on  the  angular  frequency  of  the  heat  flux  and  the  distance  from  the  surface  to 
the  first  sensor,  Fo~ 1  =  x\col  a ,  and  the  individual  curves  are  for  various  values  of  the  ratio  x2/x!.  These  two 
variables  completely  determine  the  solutions.  Considering  the  surface  temperature  first,  we  see  immediately  that  the 
solution  has  the  characteristics  of  a  low  pass  filter  but  with  a  region  of  positive  gain  in  the  range  of  angular 
frequencies  from  1  to  20.  We  also  see  that  the  magnitude  of  the  gain  increases  with  the  ratio  x2/x!.  The  phase  of  the 
solution  begins  to  lag  at  angular  frequencies  above  1  except  for  x 2/x , =4  where  there  is  a  range  from  0.4-4  where  the 
phase  angle  is  positive. 

For  the  purposes  of  comparison,  figures  2  and  3  also  contain  plots  of  the  gain  and  phase  behavior  for  the  related 
problem  of  a  simple  linear  extrapolation  from  the  internal  temperatures  to  the  surface.  The  derivations  for  this  case 
are  not  given  but  can  be  obtained  by  simply  dropping  the  terms  involving  the  second  derivatives  from  (3)  and  (4) 
and  proceeding  as  before  for  (21).  It  can  be  seen  that  requiring  the  polynomial  to  match  the  second  derivatives  at 
the  measurement  points  results  in  two  orders  of  magnitude  improvement  in  frequency  response. 

The  behavior  of  the  solution  for  heat  flux  is  given  in  figures  4  and  5  and  is  substantially  similar.  In  this  case,  the 
peak  in  the  gain  increases  with  x2/xj  but  is  significantly  lower  than  the  peak  in  the  surface  temperature  gain.  In  fact, 
for  heat  flux,  the  gain  behavior  is  quite  a  weak  function  of  x2/xi  over  the  range  from  1.5-4.  The  phase  behavior  is 
also  a  weak  function  of  x2/xi  with  all  values  showing  a  roll  off  beginning  near  an  angular  frequency  of  0.5.  Relative 
to  a  linear  extrapolation,  the  polynomial  extrapolation  increases  the  frequency  response  for  heat  flux  by  three  orders 
of  magnitude. 

In  addition  to  correctly  reproducing  the  gain  and  phase  of  a  heat  flux  signal,  the  model  must  also  be  robust  in  the 
presence  of  noise.  We  now  consider  the  effect  of  the  parameter  x2/xj  on  the  propagation  of  error  from  noise  in  the 
temperature  measurements.  Equations  (13)  and  (14)  are  expressions  for  the  uncertainties  in  surface  temperature  and 
heat  flux  when  the  smoothing  filter  approach  has  been  used.  An  alternative  approach  for  deriving  these  expressions 
is  to  start  with  (3)  and  (4)  and  assume  that  the  uncertainties  in  temperature  and  rate  of  change  of  temperature  are 
known  quantities.  This  approach  yields  succinct  expressions  that  reveal  the  role  of  x2/x!.  We  assume  the  error  in 
each  measurement  is  independent  and  the  thermocouple  positions  and  material  properties  are  error  free.  The 
sensitivity  coefficients  are  simply  the  coefficients  in  (3)  and  (4).  If  we  assume  5TX  ~  <ST, ,  STl  «  5T \  and  al  ~a2, 
the  uncertainty  of  the  surface  temperature  is, 
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(22) 


^o  = 


2  2 

X?  -x, 

( x2-x, )2 


dr 2 


2  2 

*1  x2 


(5x;  -Sxlx2+5x2)_^2 


36a2  (x2  -  Xj)2 


The  uncertainty  for  the  heat  flux  is, 


dr2  + 


(5x4  +  4x,3x,  +  4x,x2  +5x2) 
36a2 (x2  -  Xj)2 


(23) 


—  x2/x1  =2.0  — x2/x1  =3.0  — x2/x1=2.3  — Linear  x2/x1  =2.3 


Figure  2  Surface  temperature  gain  for  several  values  of  x2/xi. 
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—  x2/x1  =2  —  x2/x1=3  x2/x1  =2.3  — Linear  x2/x1  =2.3 


Angular  Frequency 


Figure  3  Surface  temperature  phase  for  several  values  of  x2lxx. 


Equations  (22)  and  (23)  show  that  uncertainty  is  a  function  of  the  positions  of  the  sensors,  and  therefore  we  can 
anticipate  that  there  will  be  optimum  locations  that  minimize  uncertainty.  To  make  quantitative  predictions  we  need 
estimates  for  ST  and  ST  and  these  were  obtained  from  experimental  data  that  will  be  presented  later.  The  specific 
values  used  here  were  ST  =1  K  and  ST  =  20  K/s.  Figure  6  shows  the  uncertainty  in  surface  temperature  based  on 
(22)  and  in  fact  it  has  a  minimum  at  x2/xl=3.6,  but  the  minimum  is  quite  broad  and  any  value  from  2.3  to  5  would 
work  as  well.  This  region  corresponds  well  to  the  region  of  maximum  accuracy  identified  previously.  However,  the 
uncertainty  increases  rapidly  for  x2/xi<2  so  this  range  should  be  avoided.  The  minimum  in  the  heat  flux  uncertainty 
(Figure  7)  based  on  (23)  is  somewhat  sharper  with  the  optimum  value  2.3.  The  heat  flux  uncertainty  has  been 
normalized  by  the  magnitude  of  the  heat  flux  so  it  can  be  seen  that  the  minimum  in  the  uncertainty  is  approximately 
0.75%  of  the  measured  value.  Again,  values  of  x2/X|<2  should  be  avoided.  Based  on  these  uncertainty  estimates 
plus  the  gain  and  phase  results  described  above,  the  optimum  value  for  x2/x  is  2. 3-2. 5. 
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0.01  0.1  1  10  100  1000 

Angular  Frequency 


Figure  4  Heat  flux  gain  for  several  values  of  x2/x! 


—  x2/x1  =2  — x2/x1=3  — x2/x1=2.3  — Linear  x2/x1  =2.3 


Angular  Frequency 


Figure  5  Heat  flux  phase  for  several  values  of  x2/xi. 
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IV.  Temperature  Dependent  Properties 

In  general,  thermal  transport  properties  are  temperature  dependent  and  the  heat  equation  is  non-linear  so 
analytical  solutions  such  as  (19)  are  not  possible.  To  estimate  the  magnitude  of  the  error  caused  by  using  (2)  we 
chose  as  a  test  case  an  idealized  version  of  the  typical  test  conditions  we  see  in  our  laboratory:  a  step  change  in  heat 
flux  from  0-4xl07  W/m2.  The  MatLab  non-linear  PDE  solver,  pdepe,  was  used  to  generate  temperature  histories  for 
the  surface  and  internal  points  and  these  were  converted  back  to  heat  flux  and  surface  temperature  using  the  model. 
Two  cases  of  x2/xj  are  shown.  For  both  cases  the  sensor  nearest  the  surface  is  located  at  a  depth  of  1.905  mm  (0.075 
inch).  The  total  thickness  of  the  slab  was  41.91  mm  (1.65  inch)  and  the  material  was  copper.  The  heat  flow  was 
assumed  to  be  one-dimensional.  The  grid  contained  131  evenly  spaced  points  and  the  temperatures  were  saved  at 
one  millisecond  intervals.  Results  using  a  261  point  grid  differed  by  approximately  0.1  K  and  did  not  change  the 
model  results. 

The  ability  of  the  model  to  reproduce  the  surface  temperature  and  heat  flux  boundary  conditions  is  shown  in 
figures  8  and  9.  For  both  cases  of  x2/x1;  after  0.1  seconds  the  polynomial  model  has  converged  to  within  0.1  K  for 
the  surface  temperature  and  within  1%  of  the  heat  flux,  but  the  magnitude  of  the  overshoot  is  significantly  smaller 
for  x2/xi=2.3  than  3  as  we  should  expect  from  the  frequency  domain  gain  plots  presented  earlier. 


Figure  6  Temperature  uncertainty  showing  location  of  minimum 
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Figure  7  Heat  flux  uncertainty 


V.  An  Application  of  the  Model 

To  illustrate  the  application  of  the  model  we  present  data  taken  in  a  sub -scale,  heat-sink,  rocket  combustion 
chamber.  Type  K  thermocouples  were  embedded  in  a  copper  block  that  formed  one  side  of  a  2.54  cm  square 
channel.  The  locations  of  the  sensors  were  those  given  above  in  the  description  of  the  temperature-dependent 
properties  calculation.  The  reactants  were  gaseous  hydrogen  and  gaseous  oxygen  at  a  mass  flow  rate  ratio  of  1:5.5. 
The  propellants  were  completely  mixed  and  reacted  at  the  measurement  location.  The  pressure  of  the  chamber  was 
42  bar  (600  psig)  and  the  mean  velocity  of  the  gas  was  approximately  250  m/s  (820  ft/s).  In  figure  10  the  upper 
figure  shows  the  temperatures  measured  within  the  block  labeled  as  “Front  Side”  for  xt  and  “Backside”  for  x2.  Also 
shown  are  results  for  surface  temperature  based  on  (3)  and  uncertainty  based  on  (13)  which  is  read  on  the  right-side 
vertical  axis.  The  reference  in  the  legend  to  TsurO  and  T3  refers  to  the  data  from  the  third  of  six  measurement 
positions.  The  data  window  was  60  milliseconds  wide  and  contained  61  measurements.  The  uncertainties  increase 
in  the  regions  where  the  rates  of  change  of  temperature  are  not  constant  as  occurs  at  the  beginning  and  end  of  the 
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test.  The  uncertainties  propagate  to  the  surface  temperature  calculation  resulting  in  the  spikes  at  1.6  and  2.2 
seconds.  The  heat  flux  results  are  shown  in  the  lower  figure.  The  maximum  value  reached  was  2.46xl07  W/m2 
(15.0  Btu/in2/s).  The  uncertainty  in  heat  flux  due  to  noise  during  the  middle  of  the  run  is  approximately  0.5%  of  the 
measured  value. 


- x2/x1=3  —  -x2/x1=2.3 


Figure  8  Effect  of  temperature  dependent  properties  on  surface  temperature  prediction.  Test 
case  is  step  change  in  heat  flux  from  0  to  4e7  W/m2  (24.5  Btu/in2/s). 


- x2/x1=3  - x2/x1=2.3 


Figure  9  Effect  of  temperature  dependent  properties  on  accuracy  of  heat  flux.  Test  case  is  a 
step  change  in  heat  flux  from  0-4e7  W/m2  (24.5  Btu/in7s).  Heat  flux  converges  to  within  1% 
after  0.1  seconds. 


15 

American  Institute  of  Aeronautics  and  Astronautics 


VI.  Conclusion 


A  method  has  been  described  for  measuring  surface  temperature  and  heat  flux  based  on  two  temperature  sensors 
embedded  in  the  wall  of  a  chamber.  The  method  does  not  require  surface  junction  thermocouples,  which  are  prone 
to  failure  and  produce  noisy  signals  in  rocket  engine  flows,  and  is  well  suited  for  studies  of  the  effects  of  surface 
features  on  heat  transfer  enhancement.  A  polynomial  extrapolation  approach  is  used  where  the  function  is 
constrained  to  match  the  temperatures  and  the  second  derivatives  with  respect  to  distance  at  the  measurement  points. 
The  method  requires  only  current  values  of  temperature  and  its  rate  of  change  and  the  boundary  and  initial 
conditions  are  arbitrary.  The  algorithm  can  be  represented  as  a  low  pass  filter  and  the  gain  and  phase  behavior  have 
been  characterized.  The  placement  of  the  sensors  affects  the  frequency  cutoff  and  the  noise  response  and  optimum 
values  for  the  relative  positions  of  the  sensors  have  been  obtained.  The  method  utilizes  time  derivatives  of 
temperature  data  which  can  be  efficiently  calculated  using  a  polynomial  smoothing  filter.  The  method  is 
computationally  efficient,  requiring  approximately  100  multiply  and  add  operations  per  measurement. 
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—  Tsurf3  K  — FrontSide  — BackSide  — T3-Noise 


- HF3  — HF3-Noise 


Figure  10  Surface  temperature  and  heat  flux  from  experimental  data  obtained  in  a 
subscale  rocket  combustion  chamber 
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